Super symmetry on the Rocks 



DESY 06-193 



Markus Ahlers 

Deutsches Elektronen-Synchrotron DESY, Notkestrafie 85, 22607 Hamburg, Germany 
E-mail: markus.ahlers@desy.de 

Abstract. In 7£-parity conserving supersymmetric (SUSY) models the lightest SUSY particle 
(LSP) is stable and a candidate for dark matter. Depending on the coupling and mass of this 
particle the life time of the next-to-lightest SUSY particle (NLSP) may be large compared to 
experimental time scales. In particular, if the NLSP is a charged particle and its decay length 
is of the order of the Earth's diameter Cherenkov telescopes might observe parallel muon-like 
tracks of NLSP pairs produced in neutrino-nucleon interactions in the Earth's interior. We have 
investigated two SUSY scenarios with a long-lived r NLSP and a gravitino LSP in view of the 
observability at the IceCube detector. 



1. Introduction 

The interactions of high energy (E > 1 PeV) cosmic particles with the Earth matter are a 
probe of physics beyond the TeV scale. It was pointed out by the authors of Ref. [1] that under 
favorable conditions a hypothetic charged particle emerging from neutrino-nucleon interactions 
in the Earth's interior could give rise to an appreciable rate of muon-like tracks in Cherenkov 
telescopes. The underlying flux of high energy neutrinos is expected from photopion processes 
of high energy protons with the cosmic microwave background (cosmogenic neutrinos) [2, 3] or 
with the ambient photon gas of cosmic accelerators [4, 5]. For the observability of these new 
particles at neutrino telescopes the low production cross section and the small branching ratio 
compared to Standard Model (SM) processes has to be balanced by a sufficiently long life time 
and a low energy loss in matter. 



2. Supersymmetry with a Long-lived Stau NLSP 

Attractive candidate scenarios are supersymmetric (SUSY) extensions of the SM with 7£-parity 
conservation. In these models the next-to-lightest super-particle (NLSP) can only decay into final 
states containing the lightest super-particle (LSP) which is stable. For a gravitino LSP with mass 
77i3/ 2 the decay length L of the NLSP in units of the Earth's diameter 2R^ is approximately [6] 

( L \ ~ ( m NLSP m3 /2 \ 2 ( -^NLSP \ (1] 
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Constraints from big bang nucleosynthesis and the cosmic microwave background yield an upper 
limit on the mass of a gravitino LSP between 10 and 100 GeV [7]. Throughout our calculations 
we will assume that the NLSP is the right-handed stau (tr), the supersymmetric scalar partner 
of the right-handed tau, and treat this particle as effectively stable corresponding to a gravitino 
mass ?7i 3/ 2 > 400 keV. 
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Figure 1. The neutrino nucleon cross section 
of the SM compared to contributions from 
neutralino and chargino exchange. 
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Figure 2. The rate of parallel stau tracks 
through a neutrino telescope like IceCube in 
terms of the neutrino energy. 



Due to 7£-parity conservation neutrino-nucleon interactions will allways produce pairs of 
super-particles, which then promptly decay into long-lived tr NLSPs. It was argued in Ref. [8] 
that the rate of single stau events is subdominant compared to the flux of muons, even for the 
most optimistic case of very light squarks with masses around 300 GeV. This is connected to 
the fact that tracks of high energy staus in a Cherenkov detector will look like low energy muons 
which have a much larger rate. Instead, a promising signal could be the observation of stau 
pairs which might show up as parallel tracks in the detector due to their large boost [']. In 
the following sections we will give a short review of the calculation of this r+r flux through a 
neutrino telescope such as IceCube [9] and comment on some phenomenological aspects. 

3. Production and Propagation of Stau Pairs 

Cosmic neutrinos are exposed to SM charged and neutral current interactions in the Earth matter 
before they might undergo a supersymmetric process at some distance I to the observer. We 
take these interaction into account by an exponential attenuation of the initial diffuse neutrino 
flux F v along the line of sight —dz = n{x)dl using the Earth's nucleon density n = p/m p shown 
in Ref. [10]. The integrated column depth of the Earth, Zq, depends on the nadir angle 6 and 
the depth of the detector center below the horizon, which we fix at 1.9 km corresponding to the 
IceCube detector [9]. For the (yet unknown) high energy cosmic neutrino flux we will adopt as 
a benchmark the Waxman-Bahcall flux of neutrinos from optically thin sources given in Ref. [4]: 
E 2 U F WB (E U ) « 6x 10~ 8 cm" 2 s" 1 sr" 1 GeV (all flavor). 

The leading order supersymmetric contribution in deep inelastic scattering (DIS) off nucleons 
consists of slepton I and squark q production via chargino \ or neutralino x° exchange 
between neutrinos and quarks. The partonic cross sections can be found in Refs. [1, 8]. 
We have analyzed two different SUSY mass spectra with a long-lived tr NLSP which we 
denote by "min m", a scenario with super-particle masses just above the experimental limits 
m x = m x o = mj = 100 GeV and = 300 GeV, and "SPS 7" with SUSY masses corresponding 
to the benchmark point SPS 7 [1 !]. The resulting neutrino-nucleon cross section is shown in 
Fig. 1. Here and in the following we have used the CTEQ6D parton distribution functions 
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Figure 3. The distribution of stau pairs 
with respect to their track separation in a 
Cherenkov detector. 
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Figure 4. The rate of stau pairs with respect 
to the nadir angle 6 compared to single muon 
events. 



fj(x,Q 2 ) [12]. 

We assume that the emerging sleptons and squarks will promptly decay into the tr NLSP 
with an average energy of (E^j/Ej = 0.5 (1.0) and {E~)q/Eq = 0.3 (0.5) for the "SPS 7" 
("min m") mass scenario, respectively. This determines the average energy loss range (zr)j^~ of 
the staus which has been studied in great detail in Refs. [13, 14, 15]. To first order, the stau 
range scales linearly with the mass and is about three orders of magnitude larger than the 
range of muons with the same energy. In our calculations we have used an approximation of 

from Ref. [13] for staus with energy Ek? > 4x 10 5 GeVx (mf/150 GeV). Below this energy 
ionization losses dominate and were used in the calculation. For a given nadir angle the maximal 
range z max of an observable stau pair is then the lesser of z®, (z^)j, and {z?)q. 

The flux of stau pairs through the detector from the decay of the emerging slepton and squark 
is then given by 



dtdAdQ dE v 



/*/*/* A 2, SUSY 

£ j dEj Jdz J ^/i^^x^x^^xe-^-^". (2) 
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Fig. 2 shows the result of Eq. (2) for the two SUSY mass scenarios "min m" and "SPS 7" 
compared to the rate of single muon events through the detector. 



4. Detection Rate at Neutrino Telescopes 

Not all of the stau pairs (Eq. (2)) might be seen as separable tracks in a Cherenkov telescope 
if they emerge from interactions too close to the detector. As an estimate we use the opening 
angle 5 between the slepton and squark in DIS to calculate the separation of stau tracks in the 
detector: x = 2Ztan(<5/2). In Fig. 3 we show the distribution of stau tracks according to their 
distance in the detector. In Tab. 1 we show the absolute rate at a cubic kilometer neutrino 
telescope such as IceCube for different detector resolutions. 

The distribution of stau pairs according to the nadir angle 6 is shown in Fig. 4. It reflects the 
exponential attenuation of the flux of neutrinos depending on the integrated column depth of 



Table 1. Rates of upward-going stau pairs for two different SUSY mass scenarios and different 
detection thresholds of the track separation x. 
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the Earth (zq) in the line of sight. The steepening of the differential rate below 33° is an effect 
of the dense inner core of the Earth. Most of the stau pairs arrive from about 85° with respect 
to nadir. The suppression at larger nadir angles is due to the small column depth of the Earth, 
Zq, compared to the energy loss range of the staus. 

On completion the IceCube detector will consist of 80 strings separated by 125 m [9], where 
each string carries 60 optical modules separated by about 17 m. Since most of the stau pairs 
will arrive from about 5° below the horizon an effective detector resolution of 50 — 100 m with 
respect to their separation seems possible. 

5. Conclusions 

Long-lived charged particles from a new physics sector may produce alternative signals at 
Cherenkov telescopes even if their production cross section and branching ratio is very small. 
Attractive candidate scenarios include supersymmetry with tr NLSPs and gravitino LSPs. For 
our choice of supersymmetric mass spectra and neutrino fluxes an observation of stau pairs at a 
rate of a few per year is possible (see Tab. 1). However, current bounds [16] on the diffuse flux 
of high energy neutrinos are about one order of magnitude above the Waxman-Bahcall flux [4] 
and an improvement of the results by a factor ~10 is possible. 
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